1. Introduction 



The influence of lattice discreteness on the properties of nonlinear systems having kink solutions was 
investigated by several authors [1]-[10]. These studies have pointed out a large variety of effects, including 
modification of kink velocity and its form and leading sometimes to the pinning of the kink on the lattice. But 
on the other side numerical investigations [9] have revealed the existence of narrow solitary waves propagating 
in certain discrete lattices without energy loss. There is also a close connection of these problems with the 
complete integrability of specific models in nonlinear cotinuum and discrete systems [11]- [13]. 

Since the pioneering paper of Krumhansl and Schrieffer [14] the role played by the elementary excita- 
' tions of kink type in the thermodynamics of one-dimensional nonlinear systems was investigated by many 

On I authors ([15]- [20] and the references therein). As concerns the influence of lattice discreteness on these 

thermodynamic properties this has received a smaller attention [20]- [23]. The class of 1-D systems mostly 
discussed is described by the Hamiltonian (in the notations of CKBT [15]) 



< 
in 



" = + ^(*«-*)' + •^}v{<i>i)] (1) 



where the nonlinearity enters only through the potential 1^ (</>), assumed to have at least two degenerate 
^ ' minima with nonvanishing curvature. It is the merit of Trullinger and Sasaki [21] to have shown that the 

. first lattice corrections are taken into account if the potential y((/>) is replaced by an effective potential 

where c? = ^ is the mean width of the static kink. In the present paper we intend to complete their results 
[ and determine also the lattice corrections to the multi-kink contributions to the free energy. 

The partition function and the free energy are usually calculated using the transfer integral operator 
, _ (TIO) method. In the thermodynamic limit the lowest eigenvalues of TIO play the most important role. 

^ ■ In the displacive/continuum limit, ^ << 1, it is possible to transform the Fredholm integral equation 

Q I representing the TIO into a differential equation. This is easily done when the two-body interaction in (1) 

O ■ has an harmonic character * [17], [19]. One obtains 



Here 7 — /31Aluq, m* — (PAcqluq)'^ , is the eigenvector of TIO corresponding to the eigenvalue e„ and 

£« = + - in - A — . (Sa) 
7 « V 7 

The left hand side of (3) contains the product ■ with A and B two noncommuting operators, and 
in order to obtain the claimed differential equation we have to put it as e^' , where the operator C can be 
written in terms of A and B using the Baker-Hausdorff formula [20], [21], [22]. This formula turns to be a 
mixed expansion in the small parameters jp and ^ — ^(m*)^^ and consequently is very convenient for 
the proposed pourposes. The first lattice correction is obtained if one retaines only the terms proportional 
with jp. Contrary the opinion of Trullinger and Sasaki [21] the same result is found if one work with the 
symmetrised or non-symmetric form of TIO. It is possible to show through a set of transformations and 



* As was shown by Guyer and Miller [17] even for more complicated interactions between nearest neigh- 
bours it is possible to write a formal expression for the partition function and the free energy using a cumulant 
expansion. 
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neglecting systematically higher powers in that the eigenvalues of TIO, with the first lattice corrections 
included, can be found from the equation [21] ** 

+ = 0. (4) 

The result is very general and valid for any kind of potential function V{(j)). 

In the next section the leading terms as an asymptotic expansion of the lowest eigenvalues of (4) will be 
found in the low temperature limit, m* » 1. The method used, briefly sketched in the Appendix, allow us 
to make a clear distinction between the contributions of phonons and of the various kink sectors. Explicit 
expressions for the lattice corrections to the kink and kink-kink contributions to the free energy will be given. 
The last section is devoted to a discussion of the results in the spirit of the CKBT [15] phenomenology of 
independent phonons and renormalised kinks. As it is known [15], [16], [19] this phenomenology is exact in 
the low temperature limit of the continuum models and the results of TruUinger and Sasaki [21] and those of 
the present paper strongly support the conjecture that it remains valid even when the first lattice corrections 
are taken into account. 



2. Lattice Corrections to the Free Energy 



Using (3a) the free energy per unit length becomes 

F = Zn (ihuoj + Acolio (5) 

where eo is the lowest eigenvalue of eq. (4). 

In the low temperature limit (3 » l,{m* >> 1) there are several ways to find approximate solutions 
of this equation, all included in so called class of the improved WKB methods [15]- [20], [23], [24] (and the 
references therein). We shall use the method developed by us in a series of papers [20], [23], [24] which has 
the advantage to make a clear distinction between the various contributions to the free energy: phonons, 
1-kink, 2- kink and so on sectors. The method is based on well known results from the theory of asymptotic 
solutions of second order differential equations depending on a large parameter [25] an is presented in the 
Appendix. As we are interested in the leading term in m* the exact equation (4) is approximated by the 
comparison equation (A. 9). With Ve/j replacing the potential V in the definition (A. 8) of the constant a 
the first term in the asymptotic expansion of the lowest eigenvalue of the isolated potential well is given by 

Eo = + O(^). (6) 



24^2 

Then the corresponding contribution to the free energy density (5) writes 

f„ = i.,„(,;^„f!) + _i_(l-jL) (7) 

and is easily identified with the first terms in the series expansion in powers of I'^/cf of the exact free energy 
of an independent gas of phonons with the dispersion relation 

oj'ik) = oJ'Jl + 4^s^n'^). (8) 



** The discrepancy between the results of [21] and [20] comes from an inadequate treatment of a certain 
equation obtained in an intermediate step of the proof in our previous paper [20] . 
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Higlica- order terms in (m*)~^/^, neglected in (6), will correspond to anharmonic interactions between the 
phonons [23]. 

The presence of the other degenerate minima determines a symmetrical splitting of each eigenvalue of 
the isolated potential well. In the case of the 0* model we obtain Eq io = Eo± to where 

^o = ^fl-^V (9) 



with c — 1. As is explained in the Appendix the value of v is determined from the boundary conditions 
(A. 11) imposed on the eigenfunctions of the comparison equation. For the model we get 

-1 = 2x^re — exp 2XJ + v{l + ln2) - -(1 + 2v)ln{l + 2v) \ . (10) 

'-{^^l \ 2 / 

This expression is valid in leading order in (m*)~^/^, but in any order in and is the starting point to find 
multi-kink contributions to the free energy. Here J is an integral defined by 

J = r VV^Fr^e dcp (11) 
Jo 

with the left turning point of the isolated potential well. As we are interested in the ground state the 
quantity e from (11) is connected to u by 



2Vrn* 

where eg is now independent on the lattice spacing I. It is seen that J depends on P/cP both through the 
integrand and the upper limit of integration. Using (12) and the expression of 14// the turning points of the 
(j)'^ model, calculated in the order 0{l'^ /(P), are given by 



Mm = 1 T \/8e~o + ^V8eo [T:^ + V^io j + O(^). (13) 

In the same order 



J = 



2^2 70 
Writing 



J = Jo + ^ Ji (15) 



where now Jo and Ji doesn't depend any more on I /<r, we get 



Jo = ^{l^l-W,-cp^) 



+ 



2V2. 

8eo 1-202 



4V2J0 v'(M?-0^)(Mi-0^) 



These integrals can be calculated in terms of complete elliptic integrals of modulus in the vicinity of unity. 
In the low temperature limit, keaping terms up to the order eo we find 

2 8 
2v^Jo ~ --?o(l+/n— ) 

3 eo 



2\/2Ji ~ -4 + 5eo. (16) 
15 



Introducing these results into (10) we get 



e-" '"(12^°') e-"^ (17) - . 



r(-i.) 



where = | Acqujo is the energy of the static kink. 

In the sine-Gordon case (periodic potential) each eigenvalue of the isolated potential well is symmetrically 
splitted in an allowed band of width 2to where to is given by (9) with c = 2, and v results from boundary 
conditions imposed on the eigenfunction at the point (j) = n. Finally we get 



-1 = 2V7re — ' exp 2XJ + 2v{l + ln2) -— ln{l + 4i/) , (10a) 

r(-zv) r(2) V 2 / 

where the integral J is defined by 

J = i VVeff-idc^. (11a) 



The integral J is calculated in the same liniar approximation in ^ and the leading terms in eo is given by 

•' = ^^((l-7S^)-|'"|-|<'-2i)) 

where eq is defined in (12) being independent on lattice corrections. Introducing (16a) into (10a) and using 
the expression E^^ = SAcqLVq for the energy of the static kink one obtains 



_!_ — «_ e >= 7^' e ^sd^ 



n-y) r(i) 

In order to find various kink contributions we have to expand the right hand side of equations (17) in 
powers of v and to write also 

V = Vk + Vkk +■■■ (18) 



The single kink term Vk comes from the expansion of 

1 



(7 = 0.577 - Euler's constant) and one obtains easily 



-u + 7Z/2 



This result is in complete agreement with that found by TruUinger and Sasaki [21]. 



In finding the kink-kink contribution v^f. we have to take into account that the expansion (18) is a 

formal expansion in powers of e~^^k ' and consequently Vkk is of the same order as v^. Then it is easily 
found that 

Vkk = vl (l + ln{l2(3Ef^) (20) - cj,^ 

One sees that Vkk gets an imaginary part satisfying the very simple relation 

Im Vkk = TTi^fe (21) 

This relation has been obtained by Zinn- Justin [26] in his analysis of the multi-instanton contributions 
in quantum mechanics and is tightly related to the non-Borel summability of the Rayleigh-Schrodinger 
perturbation expansion when the potential has degenerate minima. 
For the sine-Gordon model one obtains in a similar way 

Vkk = 2vl (7 + In 4/3^°) - ^ - |) (20) - SG 

Both (19 - 4>^) and (20 - SG) in the limit ^ ^ are in complete agreement with previous calculations of 
the kink- kink sector contributions [26]- [28]. 



3. Concluding Remarks 



According to CKBT phenomenology the thermodynamic properties of the systems described by the 
Hamiltonians of the form (1) are influenced by the existence of the static kinks [15]-[19]. A complete 
agreement between this phenomenology and the exact results of TIO in the low temperature limit is found 
if one takes into account the scattering of phonons on the static kink, leading to a renormalisation of the 
kink energy. 

The results of Trullinger and Sasaki and of the present paper are showing that lattice corrections are 

easily included into the thermodynamics of this systems. Although a complete proof of a similar CKBT phe- 
nomenology doesn't exist at present, the existing results strongly support the idea that their phenomenology 
is still valid at low temperatures. This comes both from the phonon part, which, as mentioned above, re- 

,2 

produces exactly the first terms in the scries expansion in ^ of the free energy of a phonon lattice gas, and 
from the kink contribution. As is seen from (19) in the kink contribution appear a lattice corrected kink 
energy 

= 4"' - (21) 

where e]^^ is the known unperturbed static kink energy and the correction E^}'' for the SG and model 
are given by 

^ = BAujQCQ (22) - SG 

E^°^ = ^AOwoCo (22) - 

These corrected values have been obtained also by other authors [5] [6], and as will be shown below, result 
from a very simple perturbation theory. Indeed expanding the field variable <p{x ±1) in a Taylor series the 
discrete Euler- Lagrange equation transforms into a fourth order differential equation 



Ei'' = 


1 

~6 


4°^ 


= 


1 

~10 


4°^ 
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where y = x/d and V(t> = Looking for a perturbed kink solution 



d4> 

My) = ^T{y) + ^ <^i'^(y), (24) 

where is the kink solution of the continuum limit, the linearized equation determining t/)^^^ is 

~ ^^^^^ ^ ~ ^ ^ 

The homogeneous part of (25) has appeared some years ago in the study of the kink dynamics in the 
presence of perturbations and has as solution the "translation mode" (Goldstone mode) [15], [16], [29]. Then 
eq. (25) can be solved by the method of variation of constants and the obtained results are in complete 
agreement with those already existing in literature [4]- [7], [10]. Using the same Taylor expansion for the field 
variable one can write down a corrected integral of energy. Introducing in it the new expansion of 4>k (24) 
the integration is easily performed and the result is just that presented in (22). Still unsolved remains the 
problem of calculation of renormalization of kinks due to the scattering of lattice phonons on the perturbed 
kinks, renormalization which determines the factor multiplying the exponentials in (19). 



Acknoledgement The authors would like to thank Dr. Marianne Croitoru for a fruitful colaboration 
during many years. 
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Appendix 



The lowest eigenvalues of (4) can be found, in the low temperature limit m* >> 1, using asymptotic 
methods known from the theory of second order differential equation depending on a large parameter [25]. 
The equation under investigation is of the type 

1 cP"^ 

- + V{<l>)^ = e*, (Al) 



2to' 

where the potential V{4>) has at least two degenerate minima and nonvanishing curvature at the minima. 

The problem will be solved in two distinct steps. In the first, based on Langcr's transformation, one 
looks for an uniform valid expansion of the solution near the minimum. Due to the peculiarities mentioned 
above, this transformation is chosen in such a way to give the harmonic oscillator behaviour as a leading 
term. In the second step, the existance of the other minima of l^((/>) are taken into account using symmetry 
properties of the wave functions. Each eigenvalue of an isolated potential well is slightely splitted. These 
tunneling terms are directly related to the kink contribution to the free energy. 

Let yUi and /X2 be the two turning points of an isolated potential well. One passes from the variables ^E* 
and ^ to the new ones R and x. 





cix) = / v^^y(^# (^.3) 



e-V 



- (C')2 • 

Here C' = and x is still an undefined function of </>. Then (A.l) transforms into 



(A.4) 



+ X^R = SR {A.5) 



-1 



where 5 = — d^^^ '^'^d = 2m*. The function x is subjected to the requirement to be regular and 
not to vanish in the interval of interest. In order to obtain the harmonic oscillator behaviour as the leading 
term, one takes 

((')2 = 40^(1 -x^). (A6) 

This is a function with two simple zeroes and a: = — 1 is associated with the turning point (j) = jii and 
X = +1 to = ^2- Now it is possible to determine the relation between the new variable x and the old one 
(j). Integrating (A. 6) we get 

a(7r-cos-^a; + x^/l - x'^) = ( Ve-V d(f> (A7) 
where the constant a is given by 

1 rf^2 

a = - del). (A8) 

It can be easily proved that for all the domain, even in the asymptotic region, <5 ~ ) and taking into account 
that = 2m* >> 1, the first order approximation can be found from the following comparison equation 
[25]: 

^ + Aa'X^(l-x'')R = 0, (A9) 
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whose solutions arc; the functions of the parabolic cylinder. For the model, the only convenient solution 
( decreasing exponentially when x cxd ) is Whittaker's function D^^y), where 

y = 2\/aA x and y = a\—^. (^.10) 

The eigenvalues of the isolated potential well are obtained iov v = n an integer number. In leading order in 
A one finds 

e-W « -^{\+n) + O(A^). (All) 

Now, due to existence of the other degenerate minima of the potential V{4i) each eigenvalue of the isolated 
well is splitted into two very near levels. They can be found using the boundary conditions for the wave 
function and its derivative 

— U=o = , ^-U^o. (A12) 
The solution *((/)) of (A.l) is related to the solution R{x) of (A.9) by 

^{ck) = l^r^ R{x) , 

and the explicit connection formula between the old variable and the new variable x is (asymptotic region 
of x) 

1 1 
a;2~ - / ^/\e-V\ dcj) + + ln2|a;| + 0{x-'^). (A13) 

In applying the boundary conditions (A. 11) we have to use the asymptotic expansion of Dv{y) [30] 



DM - y'^e-^ - ^ e-. (^.14) 



The value of v is found from a matching between the dominant e « and the subdominant term e * . Then 
from (A.12,A.13,A.14) one obtains 

g-iTTi/ / /-All 1 \ 

-1 = 2s/^ — -exp[2\ I VV^d(f) + 1^(1 + ln2) - -{1 + 2u)ln{l + 2u) ) . (A15) 

r(— z^) V Jo 2 / 

The sine-Gordon case can be treated in a very similar way, the solution for the comparison equation (A.9) 
being, in that case, expressed more convenient in terms of Kummer's functions having a definite parity. Each 
eigenvalue of an isolated potential well is splitted into a narrow allowed band. If the potential V{(f>) has a 
27r-periodicity the index v corresponding to the lower and upper boundary of the lowest band follows from 
similar conditions as (A. 12) calculated in the point (f) = it. Finally a very similar relation with (A. 15) is 
found, namely 

r(i + i^) 



-1 = 2s/^ 



n-h) r(i) 



exp(2X J VV^d(t> + 2i/(l + in2)- i(l + 4zy)/n(l + 4zy)^ . {A.lba) 

Here by ±0 we have denoted the turning points of the isolated potential well, centred arround the minimum 
(j) = 0. Also the definition of v is slightly modified and insted of (A. 10) we have 2^ = aX— ^. To conclude we 
have to indicate the connection relation between the symmetric splitting of each eigenvalue of the isolated 
potential well, Eq ^ io = Eo± to and the quantity i/. We get 

to = -^i' (A16) 
where c = 1 or 2 for the (p'^ or the SG model respectively. 
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